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This letter presents results for the two-loop matching coefficients for the dimension eight operators 
that contribute to Higgs production via gluon fusion. The coefficients can be used to calculate the 
first correction to the infinite top mass limit to Higgs production with large transverse momentum 
at two loops. To date such processes have been studied at two loop order only in the leading term 
in the top mass expansion. These corrections become enhanced in processes with large final state 
invariant mass, typical of multijet processes. 



One significant objective of the Large Hadron Collider 
(LHC) is to produce and study the Higgs particle. To 
fully understand this sector of the standard model, it is 
important to have precise theoretical control over the ob- 
servables associated with the Higgs in a hadronic produc- 
tion environment. The dominant mechanism for Higgs 
production is gluon fusion through a top-quark loop (for 
review see [l|). Given that the process starts at one loop, 
it could easily be enhanced by new physics. Thus an ac- 
curate prediction for the cross section and other associ- 
ated observables is of vital importance, and much work 
has been done with this in mind 0, El, 0, [f|. On the 
other hand, the fact that the process starts at one loop 
also means that radiative corrections, which are known 
to be large, are more difficult to calculate. However, the 
calculations can be greatly simplified by working in an 
effective theory where the top quark has been integrated 
out. Doing so effectively reduces the order of the calcu- 
lation by one loop at the cost of introducing errors that 
are suppressed by inverse powers of the top mass. To 
date, the focus has been on the lowest mass dimension 
six operator generated, 
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The two loop result 0, 0] for the coefficient C is given by 
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For rrih < 2mt, this leading order (in inverse powers of 
mt) contribution does an excellent job of approximating 
the inclusive Higgs production rate, with errors on the or- 
der of a few percent for a light Higgs. On the other hand, 
observables requiring a large transverse momenta for the 
Higgs, like production in association with jets or its trans- 
verse momentum spectrum, will be more susceptible to 
larger power corrections. Such processes have been cal- 
culated to two loops in the infinite top mass limit, that is, 
only including the dimension six operator Eq. ([I]) Q . 
The full m t dependence for Higgs plus jet observables is 
presently only known at one loop 



10, 11, nj, m 



To 

extend these observables to two loops is of considerable 
difficulty, however, the calculation is simplified by work- 
ing in the effective theory. The mass corrections to these 



results can be included by first calculating the match- 
ing coefficients to the set of dimension eight operators, 
and then using those results to calculate the cross section. 
The goal of this paper to present the aforementioned Wil- 
son coefficients. In a forthcoming paper the results for 
Higgs plus jet cross section will be presented. 

Finally, it is interesting to note that these dimension 
eight operators are responsible for the leading order con- 
tribution to an observable that can be used as a test for 



new physics effects. In particular, it has been shown [14 1 



that the ratio of the inclusive <7j nc to cut er 
sections 
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is approximately model independent if all the new masses 
arc sufficiently large that the effective field theory is well 
behaved. The corrections to this statement arise from 
the dimension eight operators. That is 
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where C% corresponds to some linear combinations of the 
Wilson coefficients introduced below. Thus, if 6 is mea- 
sured and found to be non-zero, then whether or not one 
can conclude there must be light new particles in the 
spectrum can only be determined once one determines 
if the contribution from the dimension eight operators is 
sufficiently small. This calculation will be taken up in a 
future paper. 



THE OPERATOR BASIS 

Below (Eqs. ([5]) to ((HJ)) is the list of all possible oper- 
ators with mass dimension eight coupling gluons to the 
Higgs, consistent with requirements of Lorentz and color 
gauge invariance. After using integration by parts to 
remove any derivatives on the Higgs field, the Bianchi 
identity was used to remove any remaining relations, thus 
giving a linearly independent basis. A minimal basis con- 



FIG. 1: Representative diagrams for the Higgs-two-gluon Ver- 
tex necessary for fixing operators with two gluon Feynman 
rules. 



FIG. 2: Representative diagrams for the Three-gluon Vertex 
necessary for fixing operator the three field strength operator. 
Taking a mass derivative with respect to the top quark gives 
the low energy limit to the Higgs-three-gluon vertex. 



sists of four operators, and is given by 
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FIG. 3: Illustration of the low energy theorem. The circles 
denote mass derivatives and the higgs momentum vanishes on 
the left hand side. 



Three of the four operators couple to two gluons and 
the Higgs, while the fourth is only involved in process 
involving three or more gluons. Its color factor only in- 
cludes the antisymmetric color structure constants. Note 
that the basis includes two operators Ob and O c that can 
be traded for operators involiving quark bilincars using 
the equations of motion. However, calculating off-shell 
will allow us to utilize the Low Energy Theorem, giving 
a computationally simpler way to calculate the matching 
coefficients, as discussed below. One can then use the 
equations of motion to simplify the basis after matching. 



METHODS 

Canonical matching involves calculating in the full and 
effective theory and then taking the difference to find the 
matching coefficients. Alternatively, one can extract the 
matching coefficient by asymptotically expanding the in- 
tegrals around hard loop momenta which arc taken to be 
of order m t fl5| and ignoring other regions which would 
cancel in the matching. Calculating in this way reduces 
the amount of work involved. 

Matching onto the basis requires asymptotically ex- 
panding double boxes in the large top mass limit. The 
Low Energy Theorem (LET) for the Higgs (see [l6| for 
an overview) allows one to reduce the complexity of the 
calculation. In its basic form, the Low Energy Theorem 
states that the amplitude for the process X — > Y + H 
can be related to the process X — > Y as 



lim Ph ^ M(X -^Y + H) 
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where ph is the four momentum of the Higgs boson, and 
m qi and are the masses and couplings of the particles 
coupling to the Higgs. Diagrammatically this is shown 
in Figure 3. 

To use the LET in the matching, one first calculates 
off shell the corrections to the Higgs-two-gluon vertex 
and match onto three of the operators (This vertex has 
been investigated in the onshcll limit for total Higgs 
production[17|), typical diagrams arc found in Figure 
1. Then one calculates the top-quark contribution to 
the three-gluon vertex in QCD. Relating this quantity 
to the Higgs-three-gluon vertex in the limit of vanishing 
Higgs four-momentum, one can fix the fourth operator's 
matching coefficient. The low energy limit itself is off- 
shell, hence the inclusion of operators that vanish by the 
equations of motion. The Low Evergy Theorem approach 
eliminates the need to calculate the 135 two-loop box dia- 
grams for the Higgs-three-gluon effective vertex. Instead 
one need only calculate 57 two- loop triangle diagrams. 

The hard contribution of the integrals is obtained by 
Taylor expanding in the external momenta. The resul ting 
expansion leaves one with a sum of bubble diagrams [15[ 
which are much simpler to evaluate. An efficient method 
to accomplish the expansion is to first reduce all inte- 
grals to scalar integrals and then performing the Taylor 
expansion following Tarasov 18|, [l9( . In this method, one 
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takes an integral of the form 
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where kj and mj are the momentum and mass associated 
with the jth propagator. For each loop momentum, we 
have introduced an auxilary vector a; and the exponential 
factor e l ^'=i fer<l! . Differentiating with respect to the 
auxiliary vectors allows one to produce any numerator in 
the loop momenta in the integral from the scalar integral. 
But before differentiating to produce the desired tensor 
integrals, one passes to the a representation, so the above 
integral has the form 
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where s are the kinematic invariants formed from the ex- 
ternal momenta and auxiliary vectors. D(a) and Q(s,a) 
are polynomials in a* and s, uniquely determined by the 
topology of the diagram. As noted above, differentiating 
the integral with respect to the auxiliary vectors gener- 
ates the desired numerator. In the a-representation the 
differentiation generates a polynomial in the external mo- 
menta whose coefficients are proportional to scalar inte- 
grals having the same form as /, but with shifted space- 
time dimension and powers of propagators. This is a 
simple consequence of the fact that Q is polynomial in 
ai and the kinematic invariants. The shift in spacetime 
dimension accounts for each derivative bringing down an 
inverse power of D(a). After differentiation, one sets the 
auxiliary vectors to zero. 

Having reduced the diagram to scalar integrals, the 
Taylor expansion can be performed similarly. The Taylor 
expansion in the a representation is equivalent to writ- 
ing out the Taylor scries in the cxp(i^j^y-) factor, and 
distributing through the a integrations over the terms. 
Again this results in further shifts in spacetime dimen- 
sions and powers of propagators. 

After the reduction to scalar integrals and Taylor ex- 
pansion, one is left with bubble integrals with propaga- 
tors of arbitrary powers, and shifted spacetime dimen- 
sions. For the case of the two-loop calculations in the 
large mass expansion, all these integrals were of the form 
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where a simple analytic result is known for all i/j and 
spacetime dimension d. To insure gauge invariance of the 
final results, all effective action vertices were computed 
in the background field gauge, and renormalizcd with the 
MS scheme at the scale 2mt- 



All diagrams were generated with FeynArts,2Q] and 
then analysized within Mathematica. 



CALCULATIONAL CHECKS 

We have checked that we reproduce the known match- 
ing results for the dimension six operator to two loops. 
This check works independently of the choice of external 
states, so we reproduce the matching to the dimension 
six operator in both processes computed. 

From the calculation we can extract the anomolous di- 
mensions of the operator basis and compare it to known 
results, thus providing another non-trivial check on the 
calculation. In the method of regions, each region devel- 
ops infrared and ultraviolet divergences, but only the sum 
over the regions contains the divergences (both UV and 
IR) of the full theory 15[ . Thus UV divergences of the 
soft regions must cancel with IR divergences of the hard 
region. With knowledge of the UV divergences of the 
effective theory (contained in the anomolous dimensions 
of the effective operators), and the one loop matching, 
one can predict the IR divergences of the hard region. 
The anomolous dimensions of the effective operators have 
been computed before [21 , 22 1 (the operators considered 
there were pure QCD operators, but have the same QCD 
renormalization properties since the higgs field is a color 
singlet). Thus it becomes a simple matter to check that 
the coefficients of the logarithms in the asymptotic ex- 
pansion are the one loop matching coefficients times the 
renormalization factor needed to substract the UV diver- 
gences of the effective theory. Thus schematically if we 
have in a calculation for the hard region (HR) 



A 2 

HR = A uv * (e-J + Log{—)) + 

M 

A 2 

Bib * (e IR + Log( — )) + Finite, 
A* 

the effective field theory {EFT) then has 



EFT = C uv *(e^+Log(^)) + 



Dir * (e IR ~ + Log{—)) + Finite. 
I- 1 



Where A is the hard scale, and p is the effective the- 
ory scale. The two arc reproducing the full theory when 
Bir = —Cuv- 

The only UV logarithms of the hard region correspond 
to the top quark mass renormalization. The two types 
of logarithms are easily distinguished by the associated 
group theory factors due to the differing representations 
of quarks and gluons. 
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RESULTS 
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The effective lagrangian resulting from integrating out 
the top to this mass order is: 



L 
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X = ^ is the yukawa coupling to the top quark. One 
takes mt in what ever scheme one renormalizes the hard 
region. Then one runs the operators to the low scale 
using a scheme that is consistent with the scheme used 
in the matching. For simplicity the coefficients are listed 
in the M S scheme. 



CONCLUSION 

We have presented the order a 2 (m^/m t ) 3 Lagrangian 
coupling the Higgs directly to gluons. This basis will 
prove useful in understanding the higher order QCD cor- 
rections to Higgs production in association with jets, 
where the range of validity of the standard effective field 
theory begins to break down due to large final state in- 
variant masses. In a forth coming paper the basis will be 
used to examine the gluon induced Higgs production with 
a large transverse momentum observable at the LHC. 
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